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ABSTRACT 

A dark-matter-only Horizon Project simulation is used to investigate the environment- and 
redshift- dependence of accretion onto both halos and subhalos. These objects grow in the 
simulation via mergers and via accretion of diffuse non-halo material, and we measure the 
combined signal from these two modes of accretion. It is found that the halo accretion rate 
varies less strongly with redshift than predicted by the Extended Press-Schechter (EPS) for- 
malism and is dominated by minor-merger and diffuse accretion events at z = 0, for all halos. 
These latter growth mechanisms may be able to drive the radio-mode feedback hypothesised 
for recent galaxy-formation models, and have both the correct accretion rate and form of 
cosmological evolution. The low redshift subhalo accretors in the simulation form a mass- 
selected subsample safely above the mass resolution limit that reside in the outer regions of 
their host, with ^ 70% beyond their host's virial radius, where they are probably not being 
significantly stripped of mass. These subhalos accrete, on average, at higher rates than halos at 
low redshift and we argue that this is due to their enhanced clustering at small scales. At clus- 
ter scales, the mass accretion rate onto halos and subhalos at low redshift is found to be only 
weakly dependent on environment and we confirm that at z ^ 2 halos accrete independently 
of their environment at all scales, as reported by other authors. By comparing our results with 
an observational study of black hole growth, we support previous suggestions that at z > 1, 
dark matter halos and their associated central black holes grew coevally, but show that by the 
present day, dark matter halos could be accreting at fractional rates that are up to a factor 3 — 4 
higher than their associated black holes. 
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1 INTRODUCTION 

In the ACDM model, structures are seeded with initial fluctu- 
ations and merge to form bound, virialized dark matter halos 
that become more massive as the universe ages. Luminous galax- 
ies form a s baryonic matter cools and condenses at hal o centres 
dWhite & Ree s 1978; Fafl & Efstathiou 1980; Blumen thaTeTai] 
1 1984 ). Dense dark halos, however, often contain embedded sub- 
halos and it has been demonstrated that low mas s subhalos can 

survive in their hosts for several billion years ( ITormenl 119971 ; 
I II ' ' — 

ITormen et al.l ll998; Moore et al. 1999). One challenge for cosmo- 
logical N-body simulations is to link dark matter halos and subha- 
los with luminous gal axies ( iBower et al.ll2006l ; IConrov et alllOOd ; 
I Vale & Ostrikeill2006 h. Understanding this relat ionship has proved 
diffic ult toiemand et al. 2004; Gap eial. 2004; lNagai & Kravtsovl 
l2005h and most explanations are provided by semi-analytic models 
dWhite & FreiiS Il99ll; ISomerville & PrimactJ [l99^; iHatton et alj 



hereafter C06). Nonetheless, a vital ingredient in explaining lumi- 
nous galaxy growth in large groups and clusters is an understanding 
of how dark matter halos and subhalos accrete mass in dense envi- 
ronments. 

The standard impleme ntation of the Extended Press-Schechter 
(hereafter EPS) formalism teond et alj[T99ll ; lLacev & Coldl 19931) 
can be used to analytically compute the average mass accretion 
onto a halo of mass Mh- iMiller et all 12006.) (hereafter M06) 
showed that: 
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l2003l; iBower et al.ll2006l ; ICattaneo et al.ll200d ; icroton et al.ll200d 
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where 5c{t) is the critical density contrast above which an object 
will collapse to form a bound structure, D{z) is the linear growth 
factor and /(AIh) is a weak function of halo mass (for alternative 
analytic expressions for hal o growth derived using E PS theory, see 
iHioteUs & P0D0I0II2OO6I and lNeistein & Dekelll2008[) . Equation Q 
can in principle be used for all redshifts and halo masses, but a 
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recent simulation study by ICohn& White! ( l2008h tested it against 
the accretion histories of massive hales at 2 = 10 and found that it 
overestimated their accretion rate. 

A simplified assumption of the EPS framework inherent in 
equation l[T) is that halos accrete at rates that do not depend on their 
environment. This restrictive assumption, however, is not a predic- 
tio n of the theory and so various authors have recently relaxed 
it. ISandvik et al. implemented a multidimensional gener- 

alization of the EPS formalism and used an ellipsoidal collapse 
model where collapse depended both on the overdensity and the 
shape of the initial density field. They found only a weak depen- 
dence between halo formation redshift and halo clustering which 
was stronger for more massive halos, in disagreement with the 
repor t ed halo assembly bias in numerical simula t ions tGao et al.l 
laOOllGao & Whitel2007l : lMaulbetsch et alj200l) . lzentneil ( l200l) 
modified the EPS formalism by using a Gaussian smoothing win- 
dow function, and Desiacques (2008) allowed the density thresh- 
old to have an environment dependence, but both authors found 
that dense large-scale environments preferentially contain halos 
that form later. We are hence lacking an EPS model that is able 
to account for halo assembly bias and predict a modified analytic 
version of equation ([TJ for the halo accretion rate. Deviations from 
the EPS accretion rate are therefore expected in the highly non- 
linear regime of cluster formation at 2 < 1, as equation ^ cannot 
account for accretion onto subhalos embedded within larger halos. 

To date, several authors have defined prescriptions for com- 
puting accretion onto halos using dark-matter-only simulations: 

• IWechsler et all ( l2002h — henceforth W02 — identified the 
mass accretion histor y (hereafter MAH) of ^ 14000 halos at 2 = 
using the ART code jKravtsov et alj|l997h in a WMAPl cosmol- 
ogy. Using their algorithm, W02 found that the accretion histories 
of their present day halos were, on average, well fitted by: 

Mh(2) = Moe-"'"^>" (3) 

where Mo is the present day mass of a halo and 0(2/) is a param- 
eter which describes its formation epoch. Ignoring the slight mass 
dependencies of q(2/(A/h)) and the /(Afn) term in equation (TJ, 
it can be seen that equation ([3} is a sensible fit for W02 to have 
chosen because in the case of an Einstein-de-Sitter (EdS) universe, 
their (^Mh^ has the same Mhz dependence as equation ([T}, differ- 
ing only in normalization ( dSc/dz = 1.686 for an EdS universe). 

• Ivan d en Bosch ( 2002) used the N-branch merger tree algo- 
rithm of ISomerville & Kolatt ( 1999) and found that a two param- 
eter fit better described the MAHs of his halos, although M06 
demonstrated that this two parameter fit becomes unphysical lo- 
cally as it pred i cts th at present day halos are not accreting mass. 
Ivan den BoschI ( 12002 *) also provided a relation for a and Zf that 
can be used in equation (|3j: 



V1.43 



(4) 



but it is more common to define 2/ as the epoch at which the present 
day halo of interest had half of its present day mass: 

In 2 



(5) 



• More recently, iMcBride et ai] ( |2009|) investigated the MAHs 
of ^ 500000 halos from the Millennium simulation with AIh > 
10^^ M© and < 2 6 and found that only ~ 25% were well 
described by equation ((Sjl. They introduced a second parameter, j3, 
and showed that a function of the form: 



prov ided a better fit to the halo MAHs. 

• iFakhouri et all ilOld) used a joint dataset from the Millen- 
nium I and II simulations and found that equation ^ held across 
five decades in mass up to 2 = 15. 

These listed accretion fits only apply when averaged across 
all environments. In order to understand accretion in dense regions 
such as clusters, one must resolve substructure and design an ac- 
cretion algorithm that can account for accretion onto halos and all 
levels of substructure. The difficulties in devising such an accretion 
algorithm are two-fold: firstly, it should define a single progenitor 
for each and every (sub)halo which accurately represents that ob- 
ject at earlier epochs, and secondly, it must conserve mass (which 
becomes harder to do when one introduces subhalos). In this study, 
outputs from a high resolution dark-matter-only N-body simulation 
have been used and a new robust method for defining accretion 
onto halos and subhalos is provided, building on previous simula- 
tion studies and moving beyond EPS theory. The primary aim is to 
investigate exactly how accretion onto halos and subhalos behaves 
as a function of redshift, mass and environment. 

One way of measuring a (sub)halo's environment is to com- 
pute the two-point correlation func tion, as this yields in formation 
on halo bias or degree of clustering.'P ercival et alj ( l2003l) used four 
ACDM simulations with differing box sizes, erg values and particle 
masses, with each simulation containing 256^ particles, to exam- 
ine four different halo merger samples at 2 = 2. They found no 
difference in clustering at this redshift between the merger samples 
of halos of a given mass. We examine the clustering of halos and 
subhalos in a higher resolution simulation and test this conclusion 
at 2 ~ 2 and at lower redshifts. 

A natural corollary is then to investigate whether the dark mat- 
ter distribution alon e has any relevance to SFR/ ^alaxy downsizing 
("Cowie et al.' 1996; Brinchmann & Ellis 2000: Baue r_et al.ll2005l: 



.Bundv et al. 2006; Faber et al. 2007; Panter et al. 20 0j|) and AGN 
down s izing dCowie et al] 12 003; Steffen et al. 2003; iBarger et af] 
I2OO5I ; iHasinger et al.l |2005| ; .Hopkins et al.i l2007i) . AGN feed- 
back provides a plausible explanation of galaxy downsizing, 
and has been su ccessfull y implemented i n semi-analytic mod- 
els dBower et al.i ,2006 ; Cattaneo et al. I I2OO& ; C06) and has been 
observed as the phenomenon responsible for t he suppression of 
star formation in ellipticals in the local universe dSchawinski et al.l 
2007, 2009). AGN downsizing is less well understood and is a two- 
fold degenerate phenomenon drive n either by low mass b lack holes 
accreting at near-Eddington rates dHeckman et al.ll2004l) or bv su- 
permassive black holes accreting at low rates ( Babic et alj2007l) . 

The structure of this paper is as follows. Section |2| describes 
the N-body simulation that was used and Section [3] explains the 
accretion algorithm. Section |4| examines accretion onto halos and 
subhalos within groups and clusters and draws comparisons with 
EPS and W02. Section [5] discusses the implications of the results 
of this paper and Section|6] the final section, lists our conclusions. 
A WMAP3 cosmology has been adopted throughout with Qm ~ 
0.24, = 0.76, Oi, = 0.042, n = 0.958, h = 0.73 and as = 
0.77. All masses are in units of Mq. 



2 THE SIMULATION 

We have analyzed outputs from one of the Horizon Project sim- 
ulationi^ which used the GAD GET- 2 code dSpringell |2005|) and 



M„{z) oc (l + 2)'^e 



(6) 



^ http://www.projet-horizon.fr 
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tracked the evolution of 512'' dark matter particles within a box 
of comoving side length 100fe~^M pc in a ACDM uni verse. 

The AdaptaHOP halo-finder ( Auber t et alj|2004l) — hereafter 
AHOP — was used to detect halos. AHOP assigns a local density 
estimate to each particle comp uted using the standard SPH kernel 
jMonaghan & Lattanziol l 19851) which weights the mass contribu- 
tions from the N closest neighbouring particles {N is usually taken 
to be 20). Halos are then resolved by imposing a density thresh- 
old criterion and by measuring local density gradients. AHOP is 
an alternative to the popular friend-of-friend (FOF) halo-finder 
jPavis etalj|l985h . which groups together particles that are spa- 
tially separated by a distance that is less than typically 20% of the 
mean inter-particle separation. Recently it has been demonstrated 
that inappropriate definitions of h alo mass can introduc e large un- 
certainties in the halo merger rate jHopkins et alJlloTol) — FOF, in 
particular, significantly overestimates the halo merger rate for halos 
that are about to merge ( Genel et al. 2009), and so we avoid using 
it. For a critical qualitative comparison between AHOP and FOF, 
see lTweed et"ZI ( l2009t) . 

In order to det ect substructure we have used the Most massive 
Sub-node Method ( iTweed et al]|2009 l) - hereafter MSM - which 
successively raises the density thresholds on the AHOP halo until 
all of its node structure has been resolved. The most massive leaf is 
then collapsed along the node tree structure to define a main halo, 
and the same process is repeated for the lower mass leaves, defining 
substructures of the main halo. For detailed descriptions of alterna- 
tive subhalo-finders, like SUBFIND, see Giocoh et al. (2010). 

The output timesteps from the Horizon simulation were sepa- 
rated by 0.01 in scale factor from 2 = 99 to the present day, but we 
restricted our analysis to halos and subhalos in the redshift range 
^ z ^ 9. The mass of each particle, Mp, was 6.8 x 10** 
and halos and subhalos with a recorded accretion value contained 
at least 40 particles. The mass of a (sub)halo used in this study cor- 
responded to the total mass, Mt, detected by the halo-finder. For 
reference, the MSM algorithm resolved 223781 objects at z = 
and ^ 20% of the se objects were subhalos. The TreeMaker code 
jTweed et al]|2009l) was then used to link together all the time out- 
puts by finding the fathers and sons of every halo and subhalo. 



3 DEVISING A HALO AND SUBHALO ACCRETION 
ALGORITHM 

This section comes in three main parts. We begin by defining the 
main branch onto a given object ("object" henceforth refers to halos 
and/or subhalos). We then provide an algorithm which identifies 
objects that take part in fake mergers. The section concludes with 
an explanation of the algorithm that was used to compute accretion 
onto bound halos and subhalos. 

3.1 A simple merger 

In Fig[T] halo i and j at timestep t2 merge to form halo k and I at 
timestep ti, where t\ > t2. In order to compute the accretion rate 
onto k, one must define a 'main father' for k and various authors 
have ado pted different prescrip tions for identifying the main father 
of a halo jSpringel et alj200ll : W02). W02, for example, define the 
main father of k as the halo that contributes the most mass to k but 
require the main father's most bound particle to be part of k if the 
main father is not at least half fc's mass. These rules force each halo 
to have a single main son and a single main father. 

There is freedom to choose the main father of k as either the 




Figure 1. Halo i and j at timestep t2 merge and form two halos, k and I 
at the later timestep ti. Halo fc's main father is j and the main branch is 
shown by the solid line. 

physically most massive father or the father that contributes the 
most mass. We have found little difference between results obtained 
from using these two definitions and so we adopt the latter defini- 
tion throughout. In Fig[T] halo fc's main father is j and the main 
branch is shown by the solid line. 

3.2 Anomalous events 

Anomalous events describe halos that spatially coincide at one 
timestep and then separate at later timesteps. These halos might 
take several timesteps to form a bound merger halo or they might 
never coincide again. One must hence be careful that their accretion 
estimator accounts for accretion onto bound objects only. 

To illustrate this point further, one would naively expect 
the mass accretion rate of halo k in Fig[T] at timestep ti to be 
(A/fc — Mj)/{ti — 12) but when this is applied to all the halos 
at timestep t\ there are a larger than expected number of negative 
accretion events (halos aren't losing mass in the hierarchical halo 
growth paradigm). Physically it is perfectly possible for mergers 
to result in mass loss along the main branch, as during a merger 
process, material is stripped from bound objects. A system of ob- 
jects undergoing a merger will, however, eventually form relaxed, 
bound objects at later times and so pinpointing the time interval 
during which mass is accreted is crucial (we do not measure mass 
loss via stripping in this work). 

3.2.1 Identifying anomalous events 

Testing to see whether an object is bound is one definitive way of 
excluding such fake events and it is common practise to sum the ki- 
netic and potential energies of each object and disregard thos e ob- 
jects whose total energy is positive jMacieiewski et al1l2009 ^. We 
combine this technique with an independent anomalous detection 
method to identify unbound objects at each redshift. 

Our prescription for identifying objects participating in 
anomalous events is as follows. The fathers of an object k at 
timestep t2 are found and if object k has two or more halo fathers 
that each donate a mass Md 20A/p, then object k is flagged as 
a possible fake merger candidate. {2QMp is chosen here rather than 
the mass resolution limit of 40Mp used in later sections, because 
20A/p is a common mass resolution limit used in other simulation 
studies and it also maximises the number of possible anomalous 
events.) The sons of k are then found and if k donates a mass 
Md 20A/p to two or more halos, then it has fragmented and 
it is identified as an anomalous event candidate. In the case of 
AHOP halos in this study, which average over their environment 
and whose substructure is not resolved, this is the sole anomalous 
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Table 1. The relative importance of unbound MSM halos and subhalos 
above the threshold mass (M > 40Afp) in the 512'^ simulation. 



Redshift 


Et > 


Anomalous (,Et ^ 0) 
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0.49 


23.1 


7.85 (84.8) 


73.5 


0.23 


23.7 


8.57 (87.3) 


73.2 


0.01 


24.1 


9.14 (89.7) 


73.4 



event criterion and the same criterion is then imposed on the next 
halo at timestep t2- 

For subhalos an additional condition is imposed. Imagine that 
two halos at timestep ts, merge to form a halo, fc, which hosts a 
subhalo at the subsequent timestep t2. Halo k and its subhalo are 
then detected as separate halos at the following timestep t\ {t\ > 
t2 > tz)- This system has transitioned over three timesteps from 
two halos, to a halo and a subhalo and back to two halos again, 
and is hence an anomalous event as no merger has taken place. The 
subhalos of a given host halo are therefore also examined if the host 
does not fragment. If a subhalo at t2 donates a mass Md ^ 2QMp 
to a halo at t\ that is a different halo to the halo son of its host, then 
it is identified as part of an anomalous event, as are its subhalos 
(if it has any) and its host. The key ideas of this anomalous event 
detection method are therefore: 

• searching for channels that receive/donate at least 20Afp 
from/to two or more different halos and 

• ensuring that the host and all associated substructures are 
flagged in the case of any one of these objects being classified as 
participating in an anomalous event. 

3.2.2 Identifying unbound objects 

Table[T]assesses the relative importance of unbound MSM objects 
above the mass threshold in the simulation {M 40A/p) for each 
of the redshifts shown in column 1 (these redshifts have been cho- 
sen because the number of subhalos increases with decreasing red- 
shift in the simulation, as clusters form). The percentages in Table 
[T] express the number of objects above the threshold mass satisfy- 
ing the condition in each column as a fraction of the total number 
of objects above the threshold mass at the redshift in question, with 
the exception of the bracketed values in column 3, which show the 
fraction of anomalous events that are unbound. 

There is a positive correlation between the independently 
identified anomalous events and unbound objects, with a large frac- 
tion of the anomalous events being unbound (henceforth unbound 
refers either to an object with total energy Et or an object 
participating in an anomalous event or both). Not all objects in col- 
umn 3 have Et ^ 0, however, and so there is a small population 
of unbound objects at each redshift that would be missed if just a 
requirement of Et were imposed on every object. 

Only bound objects above the mass threshold can have a 
recorded accretion value in this study, despite ~ 38% of all the 
objects at each of the redshifts shown in Table [T] having a mass 
below the chosen threshold limit. Bound objects below threshold, 
however, are not removed from the sample and so it is possible for 
a bound object with M < 40A/p to be a main father. We therefore 
avoid biasing the accretion events in the simulation, whilst ensuring 
that only well resolved objects have an accretion value. 

Column 4 shows the fraction of objects above the mass thresh- 
old with a recorded accretion value. A very small fraction of bound 
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Figure 2. A schematic illustrating the halosub accretion algorithm that ac- 
counts for accretion onto halos and subhalos. In this example object k, 
whose main father is object j (solid line), has been identified as the main 
son of object i (solid line). The accretion onto k using the halosub method 
is therefore (1 — fj)Mf^, where fj is the fraction of k's mass that comes 
from j. Object m is not the main son of object i and because it doesn't 
have any other fathers it is skipped. Object p's main father is q, hence the 
accretion onto p is (1 — fq)Mp. The halosub method therefore only ever 
records zero or positive accretion rates. 

objects with M ^ AOMp do not have a measured accretion rate 
because they do not satisfy some additional criteria imposed by the 
accretion algorithm, which we explain in the following section. 

3.3 The accretion algorithm 

In detecting substructure, ISpringel et al ] l l200lh required that sev- 
eral of the most bound particles of the main father were included 
in the main son — this was more robust than tracking the evolu- 
tion of the single most bound particle, which essentially performs a 
random walk across time. We have defined the main son as the son 
which receives the most mass from the object of interest, consistent 
with our main father definition. 

We shall henceforth refer to the algorithm that computes ac- 
cretion onto halos and associated substructures as the "halosub" 
method and it is illustrated in Fig[2] For object i at timestep t2 the 
main son k (solid line) is identified. Using our main son definition 
this means that most of i's mass goes to k and the remainder goes 
to m and p. The father that contributes the most mass to k is then 
found; in this example j is the main father (solid line). The mass 
accretion onto k is therefore (1 — fj)Mk where fj is the fraction 
of fc's mass that comes from object j. Object k is now flagged and 
the accretion onto the other sons of i, m and p, is considered. Since 
m is not the main son of i and m doesn't have any other fathers, an 
accretion value for m is not recorded and it is flagged as an orphan. 
If however one of the sons, p, of the object of interest does experi- 
ence mass accretion, we identify the main father, q, and record the 
mass accreted: (1 — fq)Mp. Object p would then also be flagged. 
To summarise, we list the principal features of the halosub method: 

• the measured mass accretion onto an object represents the sum 
of diffuse accretion (material not bound to any resolved structure) 
and merger-driven growth 

• mass loss events are considered to be zero accretion events: 
measured accretion signals in this study are never negative 

• all objects with a recorded accretion value are bound and have 
a mass M ^ AQMp 

• no distinction is made between halos and different levels of 
substructure 

Since we do not attempt to measure the mass lost from an ob- 
ject during a given time interval, the accretion rate in this study 
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can be thought of as an upper Umit. Note that objects which only 
lose mass and have a recorded accretion rate of zero are identified 
as systems where the bound main son of the object of interest has 
only one bound father. A flagged object means that either the accre- 
tion onto that object has already been accounted for or that object 
has been identified as an orphan. 

3.4 Limitations 

Other than finite mass and time resolutions which are shortcomings 
of any simulation, we consider the growth of halos and subhalos in 
a ACDM universe without a prescription for the gas physics. The 
dark-matter-only simulation satisfies the objective of this study, 
however: to determine whether halo and subhalo accretion is de- 
pendent on environment. The accretion algorithm excludes tidal 
stripping from the measured accretion rate but objects are stripped 
of mass in the simulation as they undergo mergers and this reduces 
their mass. 



4 RESULTS 

Throughout this section: 

(i) "object" refers to halos and/or subhalos. 

(ii) the mass of an object corresponds to the total mass, Mt, 
detected by the halo-finder. 

(iii) only bound objects above the mass threshold, M ^ 40Mp, 
can have a recorded accretion value. 

(iv) the measured mass accretion is the sum of diffuse- and 
merger-driven accretion: we have not measured mass loss. 

(v) /i = M /M denotes the specific accretion rate, with units of 
Gyr~^, onto an object of mass M. 

(vi) 5 = SMh /Mh, where Mh represents the mass of a halo. 

4.1 Accretion onto dark matter halos 

4.1.1 Comparison with EPS 

FigH] shows the average accretion rate onto the AHOP halos from 
the simulation as a function of redshift and halo mass. Halos with 
recorded accretion values are binned in mass at each redshift and 
the average accretion rate for each mass bin is computed. Averages 
of the corresponding mass bins over redshift then yield constant 
(Mh) values (W02 adopt an alternative technique, however, by 
binning the z = halos in mass and then averaging over all the 
accretion trajectories in each bin at each redshift). The solid lines 
show the accretion rates onto the AHOP halos using the halosub 
method, and the error bars indicate the la errors on the mean ac- 
cretion rate. The EPS predictions for each of the (Mh) bins, com- 
puted using equation (TJ, are shown as the dashed lines. 

Fig[3]shows that the simulation mass trajectories have a lower 
gradient across redshift than the EPS curves, which overestimate 
the accretion rate onto the lowest mass halos in the simulation at 
high redshift by a factor of ~ 2, and underestimate it by a factor 
of ~ 1.6 — 1.8 at 2; = 0. It is tempting to think that the enhanced 
accretion onto halos with respect to EPS theory at low redshift re- 
sults from the exclusion of mass loss in our measured halo accre- 
tion rate. However, EPS doesn't account for mass loss from halos 
either: halos only grow with time by construction. The offset with 
EPS should the refore be regarded as an offset in gradient and Fig[3] 
implies that the lLacev & Cold ( 1 19931) EPS formalism may only re- 
quire minor adjustment to reproduce the simulated trajectories. 
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Figure 3. The average halo accretion rate as a function of redshift and halo 
mass. The accretion values onto the AHOP halos are shown as the solid 
lines for each of the five (Mjy ) bins, with the errors corresponding to the 
Icr eiTors on the mean accretion rate. The EPS curves using equation {T) are 
shown as the dashed lines for each mass bin. 
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Figure 4. The total mass accretion rate onto the AHOP halos per comoving 
cubic Mpc as a function of redshift, halo mass and S (= SMh /Mh)- The 
mass bins correspond to the {Mh) bins in Fig|3] with the lower mass curves 
shifting to higher z. The dashed and thin solid lines show each halo mass bin 
decomposed into halos with 5 ^ 0.02 (minor-merger &; diffuse accretion) 
and (5 0.08 (major-merger & diffuse accretion) respectively. The thick 
solid lines show the mass trajectories integrated over all 5. 



4.1.2 The different accretion modes 

The mass accreted onto the AHOP halos in Fig[3]is the summed 
contribution of diffuse accretion events and minor and major- 
merger events, hence in Fig|4]we examine the relative importance 
of these accretion modes as a function of halo mass and redshift. 
At each redshift, the dimensionless quantity S (= SMh /Mh) was 
computed for each accretion event: the dashed lines and the thin 
solid lines show halos with 5 ^ 0.02 (minor-mergers & diffuse 
accretion) and S ^ 0.08 (major-mergers & diffuse accretion) re- 
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Figure 5. The total mass accretion rate onto the AHOP halos per comoving cubic Mpc as a function of halo mass and accretion mode (denoted by 5), shown 
for several redshifts. The linestyles have the same meaning as in Fig|4]except we also include halos with 0.02 < 5 < 0.08 shown by the dotted lines. 



spectively. The total mass accretion rate per comoving cubic Mpc 
for halos in a given mass bin and of a given 5 at each redshift was 
then computed. The thick solid lines show the total mass accre- 
tion rate per comoving cubic Mpc integrated over all 5. For a given 
linestyle, the lower mass curves shift to higher redshifts. 

At high redshift, all halos are found to accrete mass diffusely 
in high fractional events with the peak in activity shifting to lower 
redshifts for more massive halos. As the mass accreted onto the 
lowest mass halos via minor-mergers and diffuse accretion starts to 
plateau at low redshift, minor-merger and diffuse accretion activity 
onto the more massive halos starts to rapidly accelerate: low mass 
halos and non-halo material are being accreted onto larger struc- 
tures. By 2 = 0, the combined minor-merger and diffuse accretion 
signals dominate the growth of all halos. We further remark that the 
dashed curves have a similar cosmological evolution to the "radio- 
mode" integrated black hole accretion rate density curves found by 
C06 and Bower et al. ( 2006), but leave a more detailed discussion 
for Section|S4l 

We have tested the ability of the cut-in-delta method at distin- 
guishing between merger type by adopting the more classical pro- 
genitor mass ratio. Each progenitor j of accretor k was assumed 
to merge in turn with fc's main father i, with progenitor mass ra- 
tio X = Mi/Mj, donating fjM to accretor k at the following 
timestep, where fj denotes the fraction of k's mass that comes 



from j. Events with x ^ 3 (x > 3) were recorded as major (mi- 
nor) mergers. We found that major mergers and diffuse accretion 
events with 5 ^ 0.08 had a very similar cosmological evolution to 
the S ^ 0.08 curves in Fig|4] The minor merger and diffuse accre- 
tion events with 5 ^ 0.02 also showed a similar behaviour to the 
5 ^ 0.02 curves in Fig|4l except there were more minor mergers at 
higher redshift for all mass curves. These features do not affect our 
conclusions in Section lS!4l however. 

Fig[5] shows the shift from major-merger and diffuse- domi- 
nated growth at high redshift to minor-merger and diffuse- domi- 
nated growth at low redshift, more clearly. The linestyles have the 
same meaning as in Fig|4l except we also include the halos with 
0.02 < S < 0.08, shown by the dotted lines. It can be seen that 
minor-mergers and diffuse accretion events start to significantly 
contribute to growth for z < 0.5, and by 2 = drive accretion 
onto all halo masses. 

Qualitatively we find very similar results to Figs.|4]and[5]when 
halos are binned in jj, (= M /M) instead of 5, but the thin major- 
merger curves in Figs. |4] and |5] decouple from the thick curves at 
later epochs, for all masses. This is probably because in transition- 
ing from S to fi, one must divide S by the time interval during which 
mass is accreted, and at higher redshifts this time interval is smaller 
(time is not a linear function of redshift) and /i is hence larger than 
it is for a given 5 onto a halo of fixed mass at lower redshift. 
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Figure 6. The mean specific accretion rate onto halos and subhalos using 
the halosub method, plotted as a function of object mass for three redshifts 
corresponding to z = 0.49 (triple-dot-dashed lines), z = 0.23 (dashed 
lines) and z = 0.01 (solid lines). Lines of a given linestyle from bottom 
to top represent the accretion onto the AHOP halos and the MSM halos 
and subhalos respectively. The thick line shows the W02 result, obtained by 
using equation (3) and equation (5) ^ = 0.01. 
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Figure 7. The mean specific accretion rate as a function of mass shown for 
2 = 0.49 (triple-dot-dashed lines), z = 0.23 (dashed lines) and z = 0.01 
(solid lines) using the halosub accretion method. For a given linestyle, the 
bottom line shows the MSM halos, the middle line shows the MSM halos 
and subhalos and the top line shows the MSM subhalos. 



4.2 Accretion onto subhalos 

In this section, the AHOP halos are resolved into constituent MSM 
halos and subhalos and the halosub method is applied to these re- 
solved structures to account for accretion onto objects in groups 
and clusters. We begin by comparing the AHOP halo and MSM 
halo and subhalo specific accretion rates with the results found in 
the W02 simulation study. The mass of a halo or subhalo is hence- 
forth denoted by M, in contrast with the previous section which 
only recorded accretion onto halos with mass Mh ■ 



4.2.1 Comparing the halosub accretion algorithm with W02 

FiglH plots the average specific accretion rate for all bound ob- 
jects from the simulation as a function of average object mass 
for redshifts corresponding to 2; = 0.49 (triple-dot-dashed lines), 
z = 0.23 (dashed lines) and z — 0.01 (solid lines). These red- 
shifts have been chosen because the epoch of cluster formation is 
z < 1. The lines of a given linestyle from bottom to top represent 
the accretion onto the AHOP halos and MSM halos and subhalos 
respectively. The thick line shows the W02 result at 2: = 0.01 using 
equation ([3} (strictly, equation ^ holds at z; = but we cannot use 
our anomalous detection method at this redshift). The W02 result 
was calculated by binning in mass each z — 0.01 bound AHOP 
halo accretor and computing the corresponding average W02 a pa- 
rameter in equation (O for each mass bin (a is inversely propor- 
tional to halo formation redshift). 

The specific accretion rate onto the MSM objects is systemat- 
ically larger than the AHOP specific accretion rates at every mass 
when considering a given redshift. The MSM method resolves the 
substructure that has been averaged out in the AHOP halo, so the 
main MSM host halo and subhalos are individually less massive 
than the AHOP counterpart. The offset with MSM is probably 
caused by dividing by the larger AHOP mass, and this offset in- 
creases with increasing mass because at larger masses subhalos oc- 
cupy a larger fraction of the total AHOP mass. The mass difference 
between AHOP and the main host MSM halo therefore increases 
with increasing AHOP mass (and there are more detected halos 
than subhalos at a given redshift in the simulation, so the halos 
dominate the MSM halo and subhalo accretion signal). 

W02 fitted the accretion trajectories of their z — halos aver- 
aged over environment in a WMAPl cosmology and so their result 
can be directly tested against the AHOP curve at z — 0.01 which 
also averaged over environment, but in a universe with a WMAP3 
cosmology (W02 argue that their fitting formula does not depend 
on the chosen cosmology). We find that the W02 specific accretion 
rate has a stronger mass dependence than found for the AHOP ha- 
los in this study and so for the large galaxy- and group- sized dark 
halos, overpredicts the specific accretion rate by a factor of ~ 1.5. 

Recent studies have shown that some halo-findin g algorithms 
can lead to large uncertainti es in the halo accretion rate jOenel et aU 
I2OO9I : iHopkins et alJ[2oTol) . The disagreement across mass with 
W02 in Fig[6l however, does not result from differences in halo- 
finder: the AHOP algorithm i s very similar to the m odified bound 
density maxima technique of iBullock et alj ( l200lh used in W02. 
The disagreement most likely arises because W02 impose differ- 
ent criteria to identify the main son and main father. They adopt 
a policy, in some cases, of tracking the single most bound particle, 
which is misleading as the trajectory essentially performs a random 
walk across time. By constrast, we rigorously identify false merger 
candidates and adopt an accretion algorithm that tracks channels 
which donate/receive the most mass (and recall that by allowing 
a bound object below the mass threshold to be a main father, we 
do not bias the accretion events). Our method hence avoids using 
ad-hoc criteria. 



4.2.2 Accretion onto MSM halos and subhalos 

Fig|7] shows the specific accretion rate from bottom to top of MSM 
halos, MSM halos and subhalos, and MSM subhalos with the 
linestyles having the same meaning as in Fig[6l The average spe- 
cific accretion rates onto halos (hh) and subhalos (/is) have weak 
mass dependencies for each of the redshifts shown: (fJ-n) oc A4''''^ 
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Figure 8. The two-point correlation function plotted as a function of the mean inter-object separation (in physical coordinates) for four redshifts corresponding 
to z = 2.03,0.49,0.23 and z = 0.01. The lines in each panel represent the bound MSM halo-halo (solid), halo-subhalo (dotted) and subhalo-subhalo 
(short-dashed) accretors. The vertical long-dashed lines represent an estimate of the resolution limit in r at each redshift. 



and (/is) oc M"'^ z — 0.01, for example. Each of the halo, 
halo and subhalo, and subhalo curves shift downwards with de- 
creasing redshift: the average specific accretion rate onto a subhalo 
at z = 0.49 is a factor of 1.3 — 1.4 greater than at z = 0.01, 
for example. Major merger and diffuse accretion events at higher 
redshifts, when the universe was more dense, are more prominent. 

Fig|7]also reveals that the subhalo accretors (and this includes 
the subhalos with a zero accretion rate) accrete at a larger rate, on 
average, than the halo accretors for z < 0.5 at the mass scales 
shown. This, however, only causes a modest shift from the halo 
curve to the halo and subhalo curve at each redshift, because there 
are more halo accretors than subhalo accretors in the simulation, 
indicating that the subhalos are not responsible for the AHOP to 
MSM shift in accretion at each redshift in Fig|6] The enhanced ac- 
cretion onto subhalos can be understood by examining their mutual 
clustering and the relative velocity of their progenitors compared to 
their internal velocity, and both of these processes are discussed in 
the following sections. 



4.2.3 The clustering ofhalos and subhalos 

The main aim of this study is to investigate whether there is a re- 
lationship between the rate at which objects accrete mass and their 



environment and so in this section the clustering properties of halos 
and subhalos at different redshifts are examined. In the following 
section we specifically target accretors in different cluster-scale en- 
vironments. 

Fig[8] shows the two-point correlation function, ^, for the 
MSM accretors from the simulation as a function of the physical 
separation distance r, at the same three redshifts shown in Figs. [6] 
and [7] and at a much higher redshift of ~ 2. The lLandv & Szalavl 
1 1 19931) 11)4 estimator was used to compute ^, requiring random cat- 
alogues for each redshift. Our catalogues sampled 300000 objects 
at each redshift and were hence larger than the corresponding total 
number of detected halos and subhalos {z — 2.03 : 156120; z — 
0.49 : 211537; z = 0.23 : 216232; 2 = 0: 223781). For each 
panel in FigjS] the solid lines represent the halo-halo pairs, the 
dotted lines represent the halo-subhalo pairs and the dashed lines 
represent the subhalo-subhalo pairs. Only the clustering of bound 
accretors was measured: halo-subhalo pairs correspond to the clus- 
tering of all bound halo accretors with all bound subhalo accretors, 
for example. The vertical dashed lines show the average total diam- 
eter of an object at the redshift in question and represent an estimate 
of the resolution limit in r. 

Figd] demonstrates that subhalo-subhalo pairings are a factor 
of ~ 2 more clustered than halo-halo pairings at large physical 
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Figure 9. Distributions of Av/vc for the lialo (thick lines) and subhalo (thin lines) accretors, where Av represents the relative velocity between an accretor's 
main father and one of its other progenitors, and Vc is the accretor's circular velocity. The halo and subhalo accretors are divided into different mass bins, 
shown by the ranges of M[j and Mg respectively, and correspond to the same objects shown in Fig|2] 



scales at low redshift. This factor increases to ~ 10 — 15 at lower 
separation scales: subhalos, by definition, reside within halos and 
so cluster more strongly at small scales. The drop-off in clustering 
amplitude at the lowest scales should be ignored as this occurs at 
scales that are below the estimated resolution limit. 

The subhalo-subhalo correlation function is the sum of two 
terms: the first describes the clustering of subhalos within the same 
host and the second describes the clustering of subhalos that belong 
to different hosts. For small separations, the subhalo-subhalo cor- 
relation function has a strong contribution from pairs of subhalos 
in the same host. The clustering of halo-halo pairings is lower at 
these scales because these scales approach the size of halos, and 
so it is less common to find two halos close to each other with- 
out one or both member(s) of the pair being a subhalo. At larger 
scales, subhalos belonging to different hosts contribute strongly to 
the subhalo-subhalo clustering strength. 

The clustering amplitudes of the three curves also evolve with 
redshift: the correlation length of the subhalo-subhalo curve in- 
creases by a factor ~ 3 towards 2 = 0, for example. This is prob- 
ably because at lower redshift there are more dense clusters and 
more subhalos within a given host in the simulation, hence there is 
a stronger contribution to the subhalo-subhalo clustering amplitude 
than at higher redshift at the separation scales shown. 



4.2.4 Measuring the relative velocities between the accretors' 
progenitors 

Having established that subhalos at sub-cluster scales are more 
clustered than halos, especially at small scales, we now examine 
the distributions of Av/vc, where Av represents the relative ve- 
locity between an accretor's main father and one of its other pro- 
genitors, and Vc is the accretor's circular velocity. If Av/vc tends 
to be smaller, on average, for subhalo accretors than halo accretors 
for example, then accretion onto halos will tend to be more sup- 
pressed than accretion onto subhalos. Fig|9]shows the distributions 
of this ratio for halos (thick lines) and subhalos (thin lines) at the 
same redshifts shown in Fig|7] The Av/vc ratio was computed for 
each progenitor k (not equal to the main father j) of a given accre- 
tor: each particle accreted from the background was counted as an 
individual relative velocity event, as was each halo/subhalo progen- 
itor. So if, for example, an accretor has a main father j, a father k, 
and also accretes two particles from the background, m and n, then 
three separate relative velocities with respect to j are computed for 
that accretor. The accretors were binned in mass, and the different 
halo and subhalo mass bins are shown by the ranges of Mr and 
AIs in Fig|9] respectively. 

It can be seen from Fig|9] that the distributions of Av / Vc 
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for the halo and subhalo accretors are similar: they depend quite 
weakly on mass and their peaks coincide. 



4.2.5 Revisiting the enhanced accretion onto subhalos in F/gjT] 

It is well established that in simulations, after infall, subhalos ex- 
peri ence mass loss via tidal stripping, tidal heating and disk shock- 
ing ilGnedinet al.lll999l : lOekel et alj|2003l : iTavlor & Babuill2004 
iD'Onghia et alJIlOld) , and have a large velocity dispersion that 
scales with their host's mass. Mass stripping from an object in this 
dark-matter-only study is recorded as zero accretion, and so one 
would perhaps expect subhalos to be accreting at low rates, on av- 
erage. We have found, however, that most of the subhalo accretors 
in the simulation at 2; < 0.5 reside in the outer regions of their 
host, with ~ 70% located beyond their host's virial radius. (The 
halo virial radius roughly corresponds to r2oo, which encloses the 
region within which the halo density is at least 200 times the criti- 
cal density of the universe.) Most of these subhalos have therefore 
probably not been significantly stripped of their mass. Infact, we 
find the opposite trend in Figl?] subhalos of a given mass in the 
simulation have a larger rate of accretion, on average, than halos 
of the same mass. Having demonstrated that there is no significant 
difference between the halo and subhalo accretor distributions of 
Av/vc, we conclude that the enhanced subhalo accretion rates are 
driven by the very frequent interactions between subhalos of the 
same host at small scales (Fig[8](. Halos are less clustered at small 
scales and so accrete at lower rates, on average. 



4.3 Halo and Subhalo environment 

In this section we specifically target the effect an object's en- 
vironment at cluster scales has on the rate at which it accretes 
mass. There are two popular, independent measures of environ- 
me nt in the literature; the overdensity SrIx) in a s phere of radius 
R jLemson & Kauffmann"l999': 'Wang et al.l l2007l) and halo bias 
Isheth & Tormen 2004; Gao & White 200^^. We adopt two similar 
measures of an object's environment: the first defines an environ- 
ment mass within a cluster-sized sphere and the second uses the 
two-point correlation function. 



4. 3. 1 Environment mass 

We have defined the environment of a halo and a subhalo as the 
total mass. Me, contained within a sphere of radius R centred 
on the object of interest. Me includes the mass of all those ob- 
jects whose centres lie within the sphere as well as the mass of 
the object the sphere is centred on. We consider spheres of radii 
R = 1.46/i^^Mpc and R — 3.65/i~^Mpc because a) these scales 
represent both typical clusters and much larger clusters and b) vari- 
ous authors have found that the dependence of some halo properties 
on environment, such as halo formation redshift, are sensitive to the 
choic e of sphere radius ILemson & Kauffmann 1999; Harker e t al.l 
I2OO6I ; iHahn et alj|2009l) . Both of these environment mass defini- 
tions are applied to each bound accretor at the redshift under con- 
sideration, with only bound accretors having a recorded Me value. 
Unbound objects and resolved objects with M ^ 40A/p are not, 
however, excluded from the sample as these objects could be part 
of a bound object's environment. 

The first row of FigfTOlplots the specific accretion rate onto ha- 
los and subhalos as a function of average object mass (A/) and aver- 
age environment mass (Me) for z = 0.49 (first column), z = 0.23 



(second column) and z = 0.01 (third column) using a sphere radius 
of 1.46h~^Mpc. The second row of FigfTOlshows the results using 
a larger sphere radius of 3.65/i~^Mpc at the same three redshifts. 

The solid lines represent the environment mass bins which 
from bottom to top for the first row are: Me < 10^^'^Mq, 
IO^ '^Mq s; Me < lO^^ '^Mo and IO^^'^Mq ^ Me < 
10^^'^Mq. The triple-dot-dashed line shows the largest environ- 
ment mass bin of lO^ '^M© Me < W^^'^Mq. For the 
larger scale environments in the second row (from bottom to top): 
Me < W^'^-^Mq, IO^^-^Mq ^ Me < W^^'^Mq (solid lines) 
and lO^^-^Mo ^ Me < lO^ '^M© (triple-dot-dashed lines). The 
vertical arrow shows the direction of increasing environment for 
all panels, with the exception of the largest environment mass bins 
in the first row, which mostly lie beneath the second largest envi- 
ronment bins. The stars in each panel represent the accretion onto 
MSM halos and subhalos independent of their environment and the 
squares joined by solid lines show the EPS results. 

The relationships found in the previous sections are preserved 
in Fig[TO] the specific accretion rate increases with object mass for 
objects in most environments and decreases towards z — (as was 
shown in Fig|7]l, and EPS consistently underestimates the mass ac- 
creted onto all object masses (as was shown for halos at 2: < 1 
in FigO. The most striking feature of Fig[TOl however, is that ob- 
jects of a given mass residing in more massive environments do 
not accrete at a particularly enhanced rate compared with objects 
of the same mass in much lower mass environments. This suggests 
that the specific accretion rate onto halos and subhalos does not de- 
pend strongly on environment at cluster scales. Objects in cluster 
mass environments shown in the first row (triple-dot-dashed lines) 
mostly accrete less mass than in lower mass environments, but the 
number of objects in cluster mass surroundings is limited by the 
choice of sphere radius. This effect is not seen for the larger-scale 
environments shown in the second row, for example, where merg- 
ing between subhalos on the outskirts of the host halo is proba- 
bly driving accretion (but only at a slightly higher overall rate). 
The second row shows that the specific accretion rate only depends 
weakly on environment at larger scales that probe the outermost re- 
gions of clusters. This weak environment dependence in rows 1 and 
2 therefore seems to suggest that the increased interaction rates of 
halos in group- and cluster- mass environments are not sufficiently 
large enough to significantly overcome the large halo relative ve- 
locities, resulting in only a modest net increase in accretion. 

Halos dominate the accretion signals in Fig[TOl but we find 
the same trends at each of the chosen redshifts when just subhalos 
are plotted as a function of their mass and environment mass. There 
are two differences, however: the subhalos a) accrete at higher rates 
and b) reside only in larger mass environments. The subhalo curves 
have been omitted in Fig[lO]for clarity. 

Other authors have quantified environment by computing 
the o verdensity Sr in a sp here of radius R, rather than the 
mass ( Lemson & Kauffmann 1999': 'H arker et alJl2006l : lHahn et al.l 
[2009; Fakhouri&Ma 2009, 20IG). We therefore calculated a 
weighted environment density for each halo and subhalo accre- 
tor by using the standard S PH cubic spline window function 
jMonaghan & Lattanzid[l983) which weights the mass contribu- 
tions from objects close to the centre of the sphere more strongly 
than those further away. Fakhou ri & Ma (2009) showed that for ha- 
los more massive than IO^^Mq, the density of the object the sphere 
is centred on starts to dominate the contributions to Sr, and so 
the central object's contribution was therefore both included and 
excluded in two separate weighted environment density measures. 
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Figure 10. (fi) plotted as a function of object mass {M) and environment mass (Me) for MSM-detected objects. The lines in row one, from bottom to top, 
represent Me < lO^'^M©, lO^-^M© s£ Mb < lO^^-Sj^^^ Me < lO^^'^M© (solid lines) and IO^'-^-^Mq ^ Me < lO^^'-^Mo 

(triple-dot-dashed lines) using a sphere radius {R) of 1.46/i~^Mpc. The second row shows the results for a sphere radius of 3.65/i~^Mpc for environment 
mass bins Me < W^^-^Mq, 10^^-^Mq ^ Me < W^^-^Mq (solid lines) and lO^^'^ ^ Me < IQl^'^M© (triple-dot-dashed lines). The vertical 
arrow indicates the direction of increasing environment mass for both rows, with the exception of the highest environment mass bins in row one, which mostly 
lie beneath the second highest environment bins. The open squares joined by solid lines illustrate the EPS result using equation (T) and the open stars show 
accretion onto the MSM halos and subhalos, independent of their environment. Columns one, two and three correspond to z = 0.49, z = 0.23 and z = 0.01. 



When binned in environment density, the same weak environment 
dependence as in Fig[TO]was found in both cases. 



4.3.2 Clustering in different accretion schemes 

In this section we use the correlation function as an alternative 
means to Section |4. 3. ll of measuring an object's environment, ex- 
cept we do not restrict our analysis to just cluster scales of a few 
Mpc. We consider samples of objects with very similar masses at 
different redshifts and examine whether objects of a given mass 
which accrete at larger rates ha ve a lar ge r clus tering amplitude. 
This also tests the work by IPercival et alj l l2003h . who found that 
at 2 = 2 halos of a given mass accreting at different rates do not 
cluster differently. 

The z — 2.03 panel in Fig[TT] shows the correlation func- 
tion for all those objects whose mass satisfies 10^"-^ M@ sC M < 
10" '' Mq with < 0.35Gyr"^ (solid), 0.35Gyr"^ M < 
0.6Gyr~^ (dotted) and fi ^ 0.6Gyr~^ (dashed). The lower red- 
shift panels show the correlation function for objects whose mass 
satisfies W^^Mq s; M < IO^'^Mq with /i < O.lGyr'^ (solid). 



O.lGyr"^ s; M < 0.2Gyr"^ (dotted) and fi > 0.2Gyr"^ (dashed). 
The mass interval for z < 0.5 has been chosen because it lies be- 
low the break mass. A/*, in the mass function at these redshifts and 
so we do not bias fi. For comparison, the mass interval in the z ~ 2 
panel lies closer to M,t. The vertical dashed lines represent an es- 
timate of the resolution limit in the separation scale (same as the 
vertical dashed lines in Fig[8]l. 

At well resolved non-linear small scales for z < 0.5, objects 
with high specific accretion rates are up to a factor of ~ 3 more 
clustered than the lower accreting objects, whereas at larger lin- 
ear scales the difference in clustering between different accretors 
is much smaller. For the cluster-scale environments of the first row 
of FigjlOl corresponding to an r value of 2.92/i~^Mpc, there is 
a weak environment dependence, with objects of larger being 
slightly more clustered. Fig[TT] therefore provides further evidence 
that the mass accreted onto halos and subhalos of a given mass 
weakly depends on their environment at cluster scales. 

In contrast to the z < 0.5 behaviour, there is very little dif- 



ference in clustering between different accretors with 10^°'^M( 



M < lO^-^A/0 at z 



2 and this holds for both the linear and 
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Figure 11. The two-point correlation function plotted as a function of the mean inter-object separation (in physical coordinates) for the four redshifts corre- 
sponding to 2 = 2.03, 0.49, 0.23 and z = 0.01. The z = 2.03 panel shows the MSM objects whose masses satisfy 10^"-^ Mq ^ M < 10^''-^Mq with 
fj, < 0.35Gyr~^ (solid), 0.35Gyr~^ ^ fj, < 0.6Gyr~^ (dotted) and ^ ^ 0.6Gyr~^ (short-dashed). The lower redshift panels show the MSM objects whose 
masses satisfy lO^M© ^ M < lO^-^Af© with < ClGyr"! (solid), O-lGyr'^ s£ M < 0.2Gyr-l (dotted) and fi ^ 0.2Gyr-i (short-dashed). The 
vertical long-dashed lines represent an estimate of the resolution limit in the separation scale at each redshift. 



n on-linear scales show n. We therefore agree with the conclusions 
of lPercival etai] ( l2003h at z ~ 2 but show that they break down at 
z < 0.5, where there is a larger difference in clustering between 
high accretors and low accretors of a given mass at all scales. 



5 DISCUSSION 

5.1 Disagreement with EPS theory 

Despite its success at reproducing the dark halo mass function in 
simulations, we find that the analytic EPS calculation shows signif- 
icant departures from the halo accretion rates found in our simu- 
lation at both low and high redshift (Fig[3}. This simulation study, 
however, is n ot the first to rep ort disagreement with EPS theory at 
high redshift: ICohn & White! {2008) examined the accretion onto 
halos of mass Mh = 5 — 8 x IO^/i'^Mq at 2 = 10 and found 
that EPS overestimated the halo accretion rate by a factor ~ 1.5 
(using a lookback time of 50 Myrs). Fig[3] shows a similar be- 
haviour, with EPS overpredicting the accretion rate onto halos of 
mass Mh ^ IO^^'^Mq by a factor of ~ 2 at z = 8. One might 
expect EPS to overestimate accretion onto halos at high redshift 



because it assumes that collapse is spherica l and that the density 
barrier is fixed in height iLacev & Cole|[l993i) whereas it has been 
shown that allowing for ellipsoidal collapse and treating the critical 
density contrast for collapse as a free parameter bette r reproduces 
the N -body halo mass function ( Sheth et al. 2001 ; Shet h & TormenI 
I2OO2I) . This modification reduces the critical density contrast for 
collapse by a factor of \/0.7 (M06) which reduces /{Mh) in equa- 
tion l[T) by the same factor, causing a slight shift in the dashed 
curves in Fig|3]but otherwise having no effect on the redshift or 
mass dependence. 

The disagreement might arise because EPS theory is only ap- 
proximate: (i) it assumes spherical collapse, whereas halos in dark- 
matter-only simulations are triaxial; (ii) it contains no dynamical 
information, and so is unable, for example, to account for mass be- 
ing stripped from one halo and then being accreted onto another; 
(iii) it cannot account for accretion onto substructures; and (iv) it 
averages over halo environment. The latter restrictions are particu- 
larly problematic in the non-linear regime at z < 1, when accretion 
onto structures embedded within clusters is of interest (FigllOt. Re- 
cent attempts to incorporate an envir onment dependence into the 
EPS excursion set theory ( Maulbets ch et alj l2007l : ISandvik et al.l 
I2OO7I : IZentneJl2007l ; lDesiacauesll20oW could modify equation ([TJ 
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which might resuh in better agreement w ith our simulation results 
for z < I'm Fig llOllBenson et ^ ( l2005h highlighted further weak- 
nesses with the EPS formalism that coul d also ac count for the offset 
in Figs. [3] and [To] They showed that the lLacev & Cole (1993) EPS 
formula yields merger rates that are not symmetric under exchange 
of halo masses, and which do not predict the correct evolution of 
the Press-Schechter mass distribution, indicating that constructed 
EPS merger trees are fundamentally flawed. 

Despite these limitations, the gradients of the EPS curves in 
FigHJare only slightly steep er than the corresp onding simulation 
curves. This implies that the'Lacev & Cold <1993) EPS formalism 
may only require minor adjustment to agree more closely with the 
simulation trajectories across mass and redshift. 

5.2 The weak relationship between accretion rate and 
environment at cluster scales 

By quantifying accretion onto substructures embedded in groups 
and clusters, we have moved beyond the limited predictive power 
of the EPS formalism. Fig[7] demonstrates that subhalos accrete 
at larger rates than halos of the same mass, on average, in the 
simulation (by a factor of ~ 3 for the lowest mass subhalos at 
z = 0). At first gl ance this appears to contradict recent claims : 
lAnguloetalJ ( |2009|) and lHester & Tasitsiomil feoiol) . for example, 
have shown that subhalo-subhalo mergers are rare and that subha- 
los are severely stripped of mass, which probably means that the 
accretion rates onto their subhalos are likely to be low. The sub- 
halo accretors at low redshift in this study, however, form a sub- 
sample of subhalos that are safely above the mass resolution limit 
and that are mostly located at large distances from their host's cen- 
tre, with ~ 70% residing beyond their host's virial radius. (The 
halo virial radius approximately encloses the region within which 
the halo density is at least 200 times the critical density of the 
universe). These subhalos are probably not therefore being sig- 
nificantly stripped of mass, unlike the subhalos in recent studies. 
The mass-selected nature of our subhalo accretors and the different 
spatial distribution within their host are therefore the most likely 
causes of the apparent accretion rate discrepancy with the studies 
mentioned above. We have further shown that the subhalo accre- 
tors in this study are more clustered than the halo accretors at small 
scales (Fig[8l( and that there is no significant difference between the 
distributions of Av/vc, where An is the relative velocity between 
an accretor's main father and one of its other progenitors, and Vc 
is the accretor's circular velocity. The high subhalo accretion rates 
are therefore likely to be driven by the very frequent interactions at 
small scales with other subhalos of the same host. 

One might expect the accretion rate onto halos and subhalos to 
depend strongly on environment at larger, cluster-sized scales given 
the increased rate of interactions in dense environments, but only a 
weak dependence is found (Figs.llOlandll U. The subhalo accretors 
reside in only the most massive environments and probably accrete 
mostly locally from their nearby subhalo neighbours rather than 
their host, and so this is a possible explanation for their weak re- 
lationship between accretion rate and environment. One likely ex- 
planation for halos is that the increased interaction rates of halos 
in group- and cluster- mass environments are not sufficiently large 
enough to significantly overcome the large halo relative velocities, 
resu lting in only a modest n et increase in accretion at cluster scales. 

iFaldiouri&Mallioioh examined the environment dependence 
of accretion onto high mass halos (AIh > IO^^Mq) from the 
Millennium simulation and found a weak, negative correlation for 
galaxy-mass halos. We find a weak but positive dependence for all 



object masses in FigfTO] Our analysis, which extends theirs by ac- 
counting for accretion onto substructures, has a different expression 
for the mass accretion rate but we have found little difference in 
the results obtained from using the two expressions for bound ob- 
jects. The obvious source of t he discrepan c y is th erefore the method 
used to identify anomalies. lOenel et^ ( |2009|) have highlighted 
some fundamen t al prob lems with the 'stitching' algorithm used by 
iFakhouri & m3 ilOld , to remove anomalous events, demonstrat- 
ing that it can lead to a double counting of mergers and to a false 
counting of anomalous events as mergers. They show that this over- 
estimation of the merger rate is particularly problematic for minor- 
mergers. Predicting the effects that overestimating the merger-rate 
has on the accretion rate and how this varies as a function of envi- 
ronment is not trivial, but differences between the anomalous event 
detection methods could explain the difference in the sign of the 
trend between accretion rate and environment. 

The z — 2.03 panel in Fig[TT] reveals that at higher redshift 
when halos far outnumber subhalos in the simulation, the rate of 
accretion onto halos is independen t of environment, con firming the 
IPercival et all ^20031) result. The IPercival et J] ( |2003|) study ex- 
amined the difference in clustering at 2 = 2 between halos of 
a given mass accreting at different rates. They considered several 
mass intervals ranging from IO^^'^A/q Mh ^ IO^ ^'Mq to 
lO^'^M© ^ Mh ^ lO^^-^M© and concluded for each mass in- 
terval that halo accretion rates do not depend on environment at this 
redshift. We suggest that this apparent lack of enviro nment depen- 
dence arises because the halos in the IPercival et al.l ( 12003) study 
and to a lesser extent the halos considered in the z ~ 2 panel in 
Figim represent some of the most massive objects at 2 ~ 2 and 
hence have bias factors fe > 1 ( ISheth&Tormenll999h . These struc- 
tures are located at the highest peaks in the density field and so by 
computing the clustering amplitude of these objects one is essen- 
tially measuring the clustering pattern of the highest density peaks 
at this redshift. It is therefore unlikely that the highest mass halos 
experiencing different instantaneous accretion rates differ in their 
clustering. By contrast, the lower mass halos and subhalos in the 
z < 1 panels are less biased and so more closely track the cluster- 
ing of the underlying mass distribution. 

5.3 Comparing dark halo growth with black hole growth 

Under the assumption that, on average, black hole growth traces 
dark halo growth (so-called "pure coeval evolution"), M06 tested 
the predictions of equation (|TJ for the evolution of the integrated 
AGN luminosity density for z ^ 3. The coeval evolution model 
tests the hypothesis that the fractional mass accretion rate onto 
black holes and onto halos are equal (i.e. M /M is the same for 
both black holes and halos), and is consistent with the tight re- 
lation inferred betwee n black hole mass and g alaxy bulge mass 
dTremaine et al.ll2002l . but see lBatcheldoj[201(il for an alternative 
interpretation), and is easy to test. M06 found the predicted inte- 
grated AGN luminosity density to be in remarkable agreement with 
the bolometric AGN luminosity density measured using hard X-ray 
data. They also found that for 2: > 0.5 average black hole growth 
is well approximated by pure coeval evolution, but for z < 0.5 
the black hole luminosity density tails off more quickly than dark 
halo growth, and by 2 = is lower by a factor of ~ 2. They sug- 
gested that this slowdown in black hole accre tion could be related 
to cosmic downsizing (e.g. lBarger et alj2005l) . 

Their predictions for dark halo growth were, however, based 
on EPS theory. The simulation trajectories in Fig[3]show that EPS 
underestimates halo accretion for 2 < 1, and at 2 = is a factor 
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of ~ 1.5 — 2 lower for all halo masses. This implies that present 
day dark halos could be accreting at fractional rates that are up to 
~ 3 — 4 times higher than their associated black holes. However, for 
1 < z < 3, the simulated dark halo accretion trajectories in Fig[3] 
are reasonably well approximated by EPS. We therefore suggest 
the following scenario: for 1 < 2 < 3 black holes grow coevally 
with their dark hosts but for z < 1, the epoch of cluster formation, 
their growth significantly decouples from that of their hosts. 

It is still plausible that this decoupling is linked to the infer- 
ence that high mass black holes preferentially "turn off" at low red- 
shifts, leaving the remaining accretion activity dominated by low 
mass black holes (Heckman et al. 2004). The cause of such down- 
sizing is often assumed to be connected to the physics of the baryon 
component. Our study reinforces this assumption: if downsizing 
were a "whole halo" phenomenon it would be manifest in our dark- 
matter-only simulation, and its absence in our results confirms that 
we should seek an explanation in the baryons. 

5.4 Is halo accretion via minor-mergers and diffuse accretion 
the cause of radio-mode feedback? 

A number of authors have developed semi-analytic models of 
galaxy-formation that are tuned to re produce the galaxy lumi- 
nosity function at low redshift (e.g. I Sower et alj l200q : C06; 
iDe Lucia et alj|200^ . A key ingredient of these models is a low 
level of feedback from black hole accretion that arises in all galax- 
ies and which increases in importance towards low redshifts. The 
feedback mechanism has still not been identified: luminous, high 
accretion-rate AGN only form a small subset of the galaxy pop- 
ulation at low redshift and seem un l ikel y to provide the required 
feedback in all galaxies. ISower et"ak I ( 12006 1 required black holes to 
have relatively high accretion Eddington ratios, which may be in- 
consistent with observations: it seems that the accretion and an as- 
sociated outflow need to be hidden from view in a so-called "radio- 
mode". C06 have assumed that such a mode could be fuelled by 
Bondi accretion from the hot gas phase of their model, but the ob- 
servational evidence for such a mechanism has not been demon- 
strated either. 

The survey of IHo et al.l d 19971) revealed that a high frac- 
tion, over 40% of nearby galaxies rising to 50% — 75% of bulge 
systems, host low luminosity AGN (LLAGN), with the majority 
of LLAGN accreting at high ly sub-Edd ington rates in the range 
10"^ < Lboi/LEdd < 10~^. lHol ( [2005h argued that these are sys- 
tems where accretion occurs via a radiatively-inefficient advection- 
dominated accretion flow (ADAF). The accretion flow puffs up 
the inner disk and mater i al is a dvected towards the black hole 
( lNaravanll2003 : IHoI|20021 |2008[) . wi th outflow being channelled 
alone kinetic-energy-dominated jets jCoUin et al. 2003; Ho 2005, 
I2OOS I. This finding leads us to suggest that LLAGN, fuelled by 
low accretion rate ADAFs, may provide the radio-mode feedback. 

In our dark-matter-only study, the integrated minor-merger 
and diffuse halo accretion rate density curves in Fig|4] increase in 
importance towards the present day for all halo masses. This qual- 
itatively agrees with the cosmological evolution of the black hole 
radio-mode integrated ac cretion signal found for each of the differ- 
ent semi-analytic models jSower et al.l2006l : C06). We suggest that 
the periods when galaxy halo growth is dominated by low accretion 
rate minor-mergers and diffuse accretion events, are mirrored by 
low accretion rates onto their associated black holes, and that those 
in turn produce the LLAGN that may be the radio-mode required 
for the feedback models. 

The integrated accretion rate density onto black holes residing 



in galaxy-mass halos that are accreting diffusely and via minor- 
mergers at z = is also very similar to the integrated accretion rate 
density onto black holes residing in similar sized halos found by 
C06, who argue that radio-mode feedback is more effective in more 
massive systems. Our estimate for the total black hole accretion rate 
density tests the hypothesis that for black holes with mass Mbh 
residing in halos with mass Mh, 

y m,BH{z = 0) ~ a^§^ y M,.h[z = 0) (7) 

i i 

where a describes the non-linearity in the black hole - dark halo 
mass relation and the index i sums over all galaxy-mass dark ha- 
los and all black holes residing in these halos. Equation ([7} as- 
sumes that black hole growt h positively traces dark hal o growth, 
on average (recent claims bv lKormendv & Bendeill2oTll . however, 
argue that for bulgeless galaxies there is no such correlation be- 
tween black holes and their dark hosts, but the interpretation of this 
as meaning that t here is no such relatio n for all galaxies has been 
clearly refuted bv lVolonteri et al.ll201ll . In what follows we do not 
address the reliability of the assumption i n equation J7]l_but rather 
test its prediction for black hole growth). ^ Ferraresj j2002h found 
that Q — 1.65 and that galaxy-mass halos with Mh ~ 10^^ Mq 
have a black hole - dark halo mass ratio of ~ 10" '. According 
to Fig|4|these halos with SMh /Mh ^ 0.02 have a total accretion 
rate density of ~ 7.6 x 10''MQGyr~^Mpc~^ at 2 = 0, which when 
substituted into equation l[7) yields a total black hole accretion rate 
density of ~ lO^^'^Meyr^^Mpc^"'. This is very similar to the in- 
tegrated accretion rate density of ~ lO^^'^A/oyr^^Mpc^^ onto 
supermassive black holes at 2; = reported by C06. 

The 5 parameter (= SMh /Mh) is a free parameter in our 
model, but we have found that adopting the more classical progeni- 
tor mass ratio, x, to distinguish between merger type yields almost 
identical results to Fig|4] This provides confirmation that our 5 cuts 
are indeed capable of separating minor- and major- merger chan- 
nels. The S parameter is therefore probably no more unconstrained 
than X- 

We conclude that the low rates of accretion onto dark halos, 
driven by minor-mergers and diffuse accretion, may provide an al- 
ternative explanation to that proposed by C06 for the radio-mode 
feedback needed to reproduce the observed galaxy luminosity func- 
tion. The low redshift feedback phenomenon and its cosmologi- 
cal evolution may be driven by the cosmological evolution of halo 
minor-mergers and diffuse accretion rather than requiring accretion 
out of a hot gas phase. 



6 CONCLUSIONS 

Outputs from one of the high resolution dark-matter-only Horizon 
Project simulations have been used to investigate the environment 
and redshift dependence of accretion onto both halos and subhalos. 
We have developed a method that computes the combined merger- 
and diffuse- driven accretion onto halos and all levels of substruc- 
ture and find that: 

• Halo accretion rates vary less strongly with redshift than pre- 
dicted by the EPS formalism. This offset in gradient for each halo 
mass curve implies that perhaps minor adjustment to the EPS for- 
mula is required. 

• Comparison with an observational study of black hole growth 
leads us to suggest that dark halos at 2 = could be accreting at 
fractional rates that are up to 3 — 4 times higher than their black 
holes. 
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• Halo growth is driven by minor-mergers and diffuse accretion 
at low redshift. These latter accretion modes have both the correct 
cosmological evolution and inferred integrated black hole accretion 
rate density at a = to drive radio-mode feedback, which has been 
hypothesised in recent semi-analytic galaxy-formation models as 
the feedback required to reproduce the galaxy luminosity function 
at low redshift. Radio-mode feedback may therefore be driven by 
dark halo minor-mergers and diffuse accretion, rather than accre- 
tion of hot gas onto black holes, as has been recently argued. 

• The low redshift subhalo accretors in the simulation form a 
mass-selected subsample safely above the mass resolution limit and 
mostly reside in the outer regions of their host, with ~ 70% beyond 
their host's virial radius, and are probably not therefore being sig- 
nificantly stripped of mass. These subhalos accrete at higher rates 
than halos, on average, at low redshifts. We demonstrate that this 
is due to their enhanced mutual clustering at small scales: there 
is no significant difference between the halo and subhalo accretor 
distributions of Av/vc, where At; represents the relative velocity 
between an accretor's main father and one of its other progenitors, 
and Vc is the accretor's circular velocity. The very frequent interac- 
tions with other subhalos of the same host drive the high subhalo 
accretion rates. 

• Accretion rates onto halos and subhalos depend only weakly 
on environment at cluster scales. For halos, it appears that the in- 
creased interaction rates in group- and cluster- mass environments 
are not sufficiently large enough to significantly overcome the large 
halo relative velocities, resulting in only a modest net increase in 
accretion at cluster scales. The subhalo accretors only reside in the 
densest environments and they are likely to be accreting mostly 
from their nearby subhalo neighbours, rather than from their host. 
We further demonstrate that halos accrete independently of their 
environment at z ~ 2, as has been found by other authors, but 
show that this behaviour results from examining the clustering of 
the most massive halos with large bias factors. When less massive 
halos below M* at low redshift are considered, a weak dependence 
between accretion rate and environment at cluster scales arises. 
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